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PROBABILITY AND FREQUENCY CHARACTERISTICS OF 
SOME FLIGHT BUFFET LOADS 
By Wilber, B. Huston and T. H. Skopins Id. 


SUMMARY 


The frequency characteristics and statistical properties of the 
buffet loads measured on the unswept wing ’ and tail of a fighter airplane 
have been studied in the stall and in the shock regime. The results 
indicate that the wing loads in buffeting can be treated as the Gaussian 
response of a simple elastic system. The tail loads appear to represent 
a more complicated process . 


INTRODUCTION 


The fluctuating loads imposed on an aircraft structure in the stalled 
flight condition have been studied extensively since their potentially 
destructive character was revealed in 1930 (ref. 1). Although first 
interest was centered on "tail buffeting, " in which loads are induced on 
the tail by flow disturbances in the wake which have their origin in the 
separated flow over the stalled wing, it was soon realized that this 
induced or secondary buffeting was only part of a larger class of sepa- 
ration phenomena and that the separated flow over the wing could also 
produce appreciable loads in the airplane structure. The term "buffeting" 
is now generally used to include this primary buffeting as well as the 
secondary or induced type. As the speed of airplanes has been extended 
into the transonic and supersonic ranges, the term buffeting has also 
been applied when the separation is associated with shock formation or 
other complex flow patterns in the higher speed ranges. Buffeting is 
also used to characterize both the aerodynamic excitation and the result- 
ant structural loads. 

In attempts at quantitative treatment of buffet loads, a major diffi- 
culty is found in the random character of the fluctuations, which has led 
to uncertainty in the application of any simple numerical measure (such 
as the frequently used half of the peak-to-peak value) . This randomness 
has been recognized by many workers, but the possible applicability of 
the analytical and experimental tools developed in the study of random 
processes was first pointed out in 1951 (ref. 2). . Same experimental 
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evidence for this applicability was provided by the analysis of refer- 
ence 3 * where the peat values of wing and tail root structural shear 
measured in stalls of varying duration at various speeds- and altitudes 
were fo und to correlate well with certain parameters suggested hy con- 
sideration of the lin ear response of an aer odynami c ally damped elastic 
system to an aerodynamic excitation which is a Gaussian random process. 

The Gaussian random process , which is discussed in more detail in 
the body of the paper and in appendix A, is a special type of random 
process about which much more is known than any other type (refs, k 
to 7)- Of particular interest is the fact that the statistical proper- 
ties of such a process are simply related to its frequency characteristics 
as represented hy the power spectrum. A primary purpose of the present 
study, therefore, is to examine the character of the wing and tail loads 
in buffeting in order to determine how well buffet loads approximate a 
Gaussian random process. For this purpose a number of characteristics 
of the load time histories are evaluated and compared with what would 
be expected for a Gaussian process. In addition to the power spectrum, 
the characteristics examined are the probability distribution of the 
loads, the frequency of zero crossings, and the probability distribution 
of the peak loads. These properties were chosen because they could, be 
readily evaluated and are of interest in connection with the Btudy of 
fatigue. 

Inasmuch as buffeting was encountered during maneuvering flight, the 
time histories of measured wing and tail loads contained large-amplitude 
components at low frequencies in addition to the more rapid load fluctua- 
tions of buffeting. A brief account of the numerical methods utilized 
for the study of buffeting under these conditions is included. 


SYMBOIS 


Cjj airplane normal-force coefficient 

f frequency, cps 

fp folding frequency, cps 

f Q average number of crossings per second of the zero axis with 

positive (or negative) slope 

fp average number of positive (or negative) peaks per second 

constants of the peak probability distribution (defined in 
eqs. ( 9 ) and ( 10 )) 


^ 1*^2 
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L buffet component of root structural shear load, lb 

AL largest value of L in a specified time interval, 3b 

M Mach number 

W number in sample 

N T average number of peaks per second -which exceed a given value 

of L 

P(y) probability that a value of the random variable y "will be 
exceeded 


~ %B 
^BB 


penetration (beyond the buffet boundary). 


dynamic pressure, lb/sq. ft 
time, sec 

a specified time interval, sec 


a random variable 


standard variable, X ~ V 
* o 

angle of attack, deg 
standard deviation 

power spectral density function of buffet load, lb 2 /cps 
row matrix 


column matrix 


rectangular .matrix 


Subscripts: 


buffet boundary- 

relating to normal probability functions 
relating to peak probability functions 
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t tall 

w wing 

max maximum 

Time differentiation is denoted by dots (as y, y, etc.) 
SCOPE AMD GENERAL CHARACTERISTICS OF FLIGHT BUFFET-IOAD DATA 


The flight data analyzed in the present investigation were acquired 
in the course of the flight investigation of the wing and tail buffet 
loads reported in reference 3- The airplane, flight instrumentation, 
and data-reduction procedures are described in reference 3* Some char- 
acteristic structural, vibration frequencies of the airplane are listed 
in table I. The present study, which deals with time histories of the 
wing nnfl tail root structural shear, required a more detailed evaluation 
of the strain-gage records than was required for reference 3, which deals 
with may 1 mum values (one-half of the largest peak-to-peak fluctuation) 
encountered in a run. The measured root structural shear in buffeting 
was composed of two components: a maneuvering load (slow variations 

over a range of 20,000 pounds in the case of the wing) and a superimposed 
buffet load ( comparatively rapid fluctuations of as much as ± 3,900 pounds). 
For study of the statistical properties of the buffet loads, separation 
of the buffet component from the measured values of load was required. 

Thin separation was accomplished by means of a numerical filtering tech- 
nique which, since it may be of use in the study of other flight buffet 
measurements, is described in appendix B. 


Maneuvers Investigated 

Of the 194 r uns reported in reference 3> six were selected for the 
detailed evaluation required in the present investigation, three in the 
stall regime and three in the shock regime. Four runs were obtained 
with the basic airplane and two with the modified airplane, that is, the 
airplane with 100-pound weights added internally near the wing tips in 
order to lower the natural frequency of the wing in the fundam en tal 
bending mode from 11. 7 to 9.3 cps. The results for two of the six runs 
are reported herein, since they were found to be "typical. Run A is a 
gradual pull-up to the stall with the modified airplane at a Mach number 
of 0.46} run B, with the basic configuration, is a gradual pull-up into 
the buffet region at a Mach number of 0.74. The variation of airplane 
normal-force coefficient and Mach number in relation to the buffet bound- 
ary of the airplane is shown for these runs in figure 1. The arrows in 
figure 1 represent flight conditions for the onset and end of buffeting 
as evidenced by the accelerometer at the center of gravity. 
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Time History of Flight Condition and Buffet Intensity 

The time variations of airplane normal-force coefficient, Mach 
number, and buffet intensity for runs A and B are shown in' figure 2. 
Angle -of -attack values are also given for run A. The measure of buffet 
intensity, denoted by the symbol a(t), is the root-mean-square value 
of the buffet ing shear in a l/2-second period; that is. 



where T is taken as a period of l/2 second. This value was computed 
at successive overlapping 2/lO-second intervals and each value iB plotted 
in figure 3 at the time corresponding to the middle of the l/2-second 
period. A l/2-second period, involving the root mean square of 50 meas- 
ured values, was selected after same trial as being long enough to give 
a sample of reasonable stability in a statistical sense, yet short enough 
to reflect changes in level associated with changes in flight condition. 

Stall regime , run A .- For readier comparison- of the buffet intensity 
with airplane operating' condition in the stall, the data of figure 2(a) 
have been replotted in figures 3(a) and (b) as cr(t) against both Cjj 

and a. In these plots the square symbols are used to distinguish the 
data applicable to the stall recovery (t > 5 -secondB in figure 2(a)). 

No particular trend of wing or tail buffet load with Cjj is evident in 

figure 3, or of wing buffet load with- a, but a possible correlation of 
tail buffet load with a is indicated. In an attempt to decide on an 
objective basis whether the scatter shewn in figure 3 was larger than 
could reasonably be ascribed to chance, and thus might indicate the 
presence of same systematic factor in the data, a standard statistical 
test (Bartlett*s, ref. 8) for homogeneity was applied to the values 
of a(t). For the wing loads (fig. 2(a)) this test su gg ests that the 
scatter shown is within the limits to be expected on the basis of chance, 
but the scatter for the tail loads is greater than would be expected on 
the basis of chance alone. This finding, together with the apparent 
correlation of tail load with an^le of attack, tends to confirm the 
indications of reference 3 that "one or more additional parameters may 
exist which are important in determining tail loads but which are not 
disclosed by the present investigation. 

Shock regime, run B . - In figure 2(b) the time histories of cr show 
marked similarities to the time history of Cjj. These values of buffet 
intensity are compared directly with Cjj in figure ^(a). Both wing and 
tail loads increase Bharply with values of Cjj above about 0.k6. The 
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trend appears to be well represented by the least-squares straight lines 
shown which fit the wing and tail loads with standard errors of estimate 
of 3 6 and 18 pounds, respectively. In fitting these lines, the data 
represented by flagged symbols were emitted. These data ‘apply to the 
rapid return to level flight (t > 4 seconds) in run B, where conditions 
were c h ang in g too rapidly to be represented adequately by l/2-second 
averages . 

Because of extraneous factors, the values of a in figure 4(a) do 
not go to zero. For the wing data, the random errors associated with 
record reading, combined with strain-gage sensitivity (l6,000 lb/in.) , 
had a standard deviation of about 60 pounds. For the tail loads, reading 
accuracy was not a factor, but a periodic ccmponent is always present 
in the tail-load record at the propeller blade passage frequency, 86 cps. 
This component of the tail load in level flight in run B has a root-mean- 
square value of about 25 pounds, and an examination of the original 
records shows that it increases somewhat with lift. These effects tend 
to obscure the detail in resolving low levels of buffeting, but they are 
judged not to affect appreciably the larger values of cr(t) . 

Although it appears that a plot of the variation of buffet intensity 
with Cjj at a given Mach number could be used to determine the value 
of Ojj at which buffeting starts, and thus to find a point on the buffet 
boundary, application of this concept is complicated in the present 
instance by both the decrease in Mach number during the run and by the 
extraneous effects mentioned previously which have obscured the buffet 
loads of low magnitude. It does appear, however, that buffeting is pres- 
ent at values of Cjj greater than 0.46, a value which compares favorably 
with the condition for onset of buffeting indicated in figure 1. 

In order to correlate the variations of buffet load with both Cjj 
and Mach number in the shock regime, a penetration parameter P was 
defined in reference 3- A given value of Cjj at a given Mach number. 
Which penetrated the buffet region by an increment Cjj - Cj ^ , was 

expressed as a ratio in terms of the maximum possible penetration 


P = 


C K ~ Sb 

° 1 W " °%B 


( 2 ) 


This formula takes care of variations of both Cjj- and Mach number. The 

buffet loads of run B are plotted against P in figure 4(b). The varia- 
tion with P appears to be linear. Least-squares straight lines through 



H 


NACA TN 3733 


7 


the origin fit the wihg and tall load data with standard errors of esti- 
mate of 32 pounds and lk pounds, respectively. As before, I'-ta repre- 
sented by the flagged symbols were omitted in fitting the least-squares 
straight lines . 


ANALYSIS AMD DISCUSSION 
Method of Analysis 


As discussed in appendix A, assessment of the Gaussian or non- 
Gauss ian character of the measured flight buffet loads in terms of the 
usual formal definition of Gaussian random processes would involve a 
set of related probability distributions, most of which would be diffi- 
cult to determine or to visualize. In the present paper, only the sim- 
plest of these distributions is used directly for comparison with the 
measured loads. The other properties to be examined follow from the 
nature of these probability distributions, but are det ermin ed from the 
power spectrum of the loads. The power spectrum ©(f) may be regarded 
as a frequency analysis of the mean-square value of the load fluctuations. 
The equations for the properties used specifically in the present study 
follow. 

The standard deviation a or root-mean-square value is determined 
by the area under the power spectrum as 

cr 2 = T ©(f) df (3) 

J 0 

The number of crossings per second of the zero axis with positive 
(or negative) slope is an average frequency f Q where 



(k) 


The total nu mber of positive (or negative) peaks per second is f p , 

where 




( 5 ) 



8 


HACA TN 3733 


The probability distributions of the loads and of the load peaks 
are most conv enient ly expressed in terms of the so-called st an dard 

variable z = ^ which expresses the size of a given value of load L 

in relation to the st and ard deviation a. The probability that L lies 

W n (z) dL 

between L and L + dL is -=— , where 


»,(*) 


_1_ -z 2 /2 


( 6 ) 


is the normal frequency distribution of statistics or the "normal curve 
of error." The probability of exceeding, that is, Pjj(l), is expressed 

by an integral of equation (6) as the normal probability distribution 
function 


j _00 

W L (z) dz 

z=L/cr 


(7) 


a relation which may be evaluated by use of standard statistical tables. 

The probability that a peak load will exceed a given value L, that 
is, Pp(L), is expressed by a peak probability distribution function 

P p (z) where 

P p(z)= P H ^) + |e- 2/2 [l-P H (4)] (8) 

ftnii the constants k-^ and kg are functions of the ratio fo/fp > 

*i - )fi - (f 0 / f p) 2 




f o/ f P 


( 10 ) 


Equation (8) is illustrated in figure 5(a) for various values of 
the frequency ratio fo/fp* For the limiting case fo/fp = Peak 

probability distribution Pp(z) reduces to the normal distribution 
(eq,. ( 7 )). This figure has been prepared on normal-probability paper - 
paper on which the scales are so adjusted that the normal probability 
distribution gives a straight-line plot and which thus serves to illus- 
trate the relationship between the normal and peak probability 
distributions . 
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In a Gaussian random process, the first member of the aforementioned 
set of related probability distributions can be represented by the normal 
probability distribution, equation (7). The Gaussian character of the 
buffet loads in the present paper is, therefore, first assessed by com- 
paring equation (7) with the probability distribution of the measured 
load time histories taken from the flight records. In addition, the 
distribution of the peaks is compared with equation (8), and the values 
of fQ and fp determined by counting are compared with the values 
obtained from the power spectrum by using equations (4) and (5)* Inas- 
much as the power spectrum is so closely involved in these comparisons, 
the frequency characteristics of the loads measured in the present study 
are discussed prior to the probability characteristics. 


Power Spectrum of Buffet Loads 

The frequency content of the wing and tall loads measured in runs A 
and B is shown in figure 6. These power spectra were computed from the 
load time histories by use of the numerical techniques devised by Tukey, 
as out lin ed in reference 7* Since the spectra are based on sample time 
histories of finite length read at discrete time intervals 
(At =0.01 second), they represent certain necessary compromises with 
regard to frequency range, resolution, and precision. These compromises 
are discussed in detail in reference 7. The spectra, therefore, should 
be considered as estimates of the power spectral density functions of 
the buffet loads rather than true power spectral densities. 

The power spectra of figure 6 are essentially characteristic of the 
random response of a lightly damped s ingle -degree-of -freedom system. 

The location of the large peak at low frequencies in the wing-load 
spectra correlates well with the frequency of first symmetric bending, 
9-3 cps (table I) for the modified airplane (fig. 6(a)) and 11. 7 cps 
for the basic airplane (fig. 6(c)), and thus reflects the change in fre- 
quency characteristics associated with the added wing-tip weights. The 
first assymmetric bending mode and first torsion mode also appear in 
the wing spectra but these contributions are small. The larger low- 
frequency peak^in the tall spectra appears to reflect the fuselage in 
a torsion, or "tail rocking, 1 mode at 9-8 cpsj identification of the 
other low peaks in terms of known structural modes is less certain. 


Probability Distribution of Buffet Loads 

As a first step in assessing the Gaussian character of the random 
loads represented by the power spectra of figure 6, the number of values 
from the load time history that fall in various class intervals are 
given in table II. Similar frequency-distribution data are given in 
table III for the load ratio L/P , a simple transformation of the buffet 


10 


NACA IN 3753 


load s ug gested by the strong linear dependence of load on penetration 
beyond tbe buffet boundary shown in figure Mb). Ihe "loads" co lumn s 
of tables II and HI were obtained by sorting into tbe appropriate 
classes tbe p unche d cards containing tbe time histories of load or of 
transformed load. Also given in tbe tables is tbe mean value of each 
frequency dis tribution and tbe st and a r d deviation ff. 

\ 

Tbe probability of exceeding a given value of L or of L/P, based 
on the load frequency data of tables H and III* is s^srwn in figures 7 
n-pfl 8 . These figures have been prepared on normal-probability paper. 

Since the normal probability distribution Pjj(L) (eq. (7)) gives a 
straight-line plot, use of probability plotting paper constitutes one 
of the simplest, qualitative tests for a normal distribution. For the 
shock regime, run B, the distributions of both and Lfc (fig. 7) 

have exc essive occurrences of very large and very sm a l l, values as c can- 
pared with a no rma l distribution. For the stall regime, run A, the 
normal dis tribution is a representation of the data which would be ade- 
quate for most engineering purposes , although it is somewhat better for 
the wing loads than for the tail loads. The transformed loads for wing 
and tail L/P (fig. 8 ) correspond even more closely to the normal distri- 
bution t han do the stall loads of run A. 

A standard statistical test of the significance of such differences 
between an observed distribution and an assumed distribution as are 

evident in figures 7 and 8 is provided by the X 2 (chi-square) test. 

(See, for example, ref. 8 .) The numerical quantities required for the 
use of this test are given at the bottom of .tables II and III. Included 

are the value of X 2 and for comparison the value X 2 _, which is the 

.up 

value of X 2 at the 5-percent level of significance. In accordance 

with usual practice, a calculated value of X 2 less than that for the 
5 -percent level of significance may be regarded as indicating that the 
differences between the observed distribution and the normal distribution 
are not significant. On the basis of the qualitative comparison of fig- 
ures 7 8 and on the more objective basis of the chi-square test, it 

is concluded that the buffet loads in the stall regime and the transformed 
loads in the shock regime display one important characteristic of a 
Gaussian r andean process, normality of the first probability distribution. 


Probability Distribution of Peak Buffet Loads 

For comparison with the peaks of a Gaussian random process, the 
frequency distributions of the buffet-load peaks of runs A and B are 
given in table II and the peaks of the transformed loads are given in 
table TTT - The frequencies of both positive and negative peaks are 
tabulated. The distinction between positive and negative peaks is based 
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not on the relative position of the zero axis, hut on the usual defini- 
tions of the maxi mums and mln-tnuims of a continuous f unc tion - - 'I'hiiH a 
peak occurs -when the first derivative of a function is zero; it is a 
positive peak when the second derivative is negative (positive curvature) 
and it is a negative peak when the second derivative is positive. 

The peak probability distributions corresponding to the peak data 
of tables IX and 111 are plotted in figure 9* Thp abscissa in this 
figure is the standard variable z, equal to L/cr, or, in the case of 

the transformed load, The probability Pp(z) is the probability 

that a peak will exceed a given value of z. For the purposes of fig- 
ure 9 , the negative peak distributions of ^ach load were inverted about 
the zero axis and combined with the corresp onding positive peak distri- 
butions to give an average or effective peak distribution. This procedure 
was adopted to simplify the presentation and to increase the size of the 
comparatively small samples, since there appeared to be no special dif fer - 
ences between the positive and negative distributions. 

For comparison of an observed peak distribution with the peak distri- 
bution of a Gaussian random process, a value of the ratio fo/fp is 
required for use In equation ( 8 ). The lines plotted in figure 9 repre- 
sent the peak probability distribution function Pp(z) (eq. ( 8 )) for two 

estimates of the ratio fo/fp* The solid lines are for the valueB obtained 

by counting zero crossings and peaks in the time history; the dashed Hupr 
represent values obtained from the respective power spectra of the load 
by use of equations (4) and ( 5 ) . The numerical values involved are summa- 
rized in table IV. The spectrum values of f Q /fp average 0.12 nnwiler 

than the corresponding time-history values but the effect of this differ- 
ence on the shape of the peak distribution is small, especially at the 
level of the larger, more infrequent peaks. For the loads which hnfl a 
Gaussian first probability distribution (run. A and the transformed loads 
of run B) the observed peak probability distributions appear to be in 
good agreement with the theoretical distribution. For run B, which had 
loads that were not no rm al l y distributed, the peak distribution appears 
to reflect the excess of very large and very small values. Thp qualita- 
tive results of this comparison are confirmed by the chi-square test, 
for which the pertinent numerical data are given at the bottom of the 
peak frequency columns of tables II and HI. The values of X 2 shown 
for the peak loads of run A and the transformed peak loads of run B are 
not significant at the 5-percent level. These distributions are, there- 
fore, judged to display another characteristic of a Gaussian random pro- 
cess; that .is, the probability that a peak will exceed a given value iB 
given by the peak probability distribution for a Gaussian process. Thp 
test was not applied to the peak loads of run B in table II in view of 
the nonstationary character of the loads prior to transformation and the 
large departures of the peak distributions shown in figure 9 . 
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In table IV the values of and of fp ' obtained by use of equa- 
tions (4) and ( 5 ) from tbe power-spectrum estimates of runs A and B are 
compared with the values actually observed. Ihe spectrum value of the 
zero-crossing frequency exceeded the observed time-history value by 1.8 
to -0.6 cps. The estimates of peak frequency fp from the power spectra 
exceeded the time-history values by amounts ranging from 5*3 to 8.6 cps, 
or as much as 33 percent. Values of fQ are thus obtained from the 
spectrum more accurately than are the values of fp. For the purpose of 
load prediction, the frequency of the larger loads is primarily determined 
by the value of fQ, and the agreement between the values of fQ shown 
in table IV is considered satisfactory. The fact that the spectrum val- 
ues of fp ‘ are consistently larger than the observed values may perhaps 
indicate same non- Gaussian, character of the buffet process which is 
appearing in the second derivative. On the other hand. Ina ccuracies in 
the power estimates at the higher frequencies which do not greatly affect 
the precision of the value of a can be magnified by differentiation 
(eq.. (A24-)) and can thus affect the accuracy with which the spectrum of the 
second derivative can be obtained by this method. Although the apparent 
deviation from the Gaussian distribution which is evident in the second 
derivative is of considerable theoretical interest and appears to warrant 
further investigation, for most practical purposes buffeting can be 
treated as a Gaussian random process. 


Application of Present Results to Buffet Research 

Two results of the present study are of special interest in connec- 
tion with the investigation of buffeting. These results have to do with 
the Gaussian character of the loads and with the simplicity of the power 
spectral densities of the loads. 

The Gaussian character of the loads is significant in view of the 
wealth of statistical information available about Gaussian distributions 
as outlined in the section entitled "Method of Analysis" and in appen- 
dix A. As applied to buffeting, such statistical properties as the 
percentage of time a given load is exceeded, the number of cycles, and 
the magnitude of peaks are perhaps of most interest since they have 
immediate application in connection with estimation of limit loads and 
fatigue effects. 

The power spectral densities of the loads measured in the present 
study may be considered as the resultant of a random fluctuating input 
acting on an elastic system. As outlined in references 2 and 3, the 
characteristics of output, input, and system may be expressed in terms 
of three related frequency functions as 

®o(f) - A 2 (f)®i(f) 



NACA TN 3733 


13 


where 

®j(f) frequency content of input force 

® Q (f) frequency content of the output or response (structural load, 
displacement, velocity, acceleration, etc.) 

A 2 (f) square of the admittance of the system to a sinusoidal force 

Buffeting is thus expressed in terms of three pert inent quantities, 
each of "which may he investigated Ind ep endent ly. Since an airplane, 
vfbrationwise, is a very complex system, the admittance could conceivably 
he a complex expression involving many degrees of freedom. It is in thi s 
connection that the simple power spectra shown in figure 6 are significant. 
Each spectrum reflects predominantly the response of a single degree of 
freedom - in the case of the wing, the first symmetric tending mode. As 
illustrated in references 2 and 3, restriction to a single degree of 
freedom can lead to some simple analytical results. In the case of the 
present unswept wing, it appears that such a restriction, as exemplified 
in appendix B of reference 3, should permit the determination of first- 
order effects at least. In view of the bending-torsion coupling which 
is characteristic of swept wings, the admittance of swept wings is possi- 
bly more complex than that of unswept wings, and a study of the power 
spectral densities of other plan forms would be of interest. 


Correlation With Previous Results 

Ina s mu ch as the present study stems from the study of peals: buffet 
loads reported in reference 3, and was undertaken partly to investigate 
some of the assumptions amide in that analysis, a comparison of some of 
the present results with those of reference 3 is indicated. In the 
analysis of the peak wing buffet loads reported in reference 3, use was 
made of an asymptotic expression for the peak probability distribution 
of a Gaussian process, obtained in reference 4 (eq.. 3.6-11) and given 
in appendix A of the present paper as equation (A35). This expression, 
for the number of positive (or negative) peaks per second which 
exceed a value of L, is 


N l « f 0 e ~ L2 / 2a2 (U) 

Equation (ll) applies almost exactly to peaks larger than 1 st andar d 
deviation. Alternatively, the peak value AL, which will occur only 
once in a period T, is obtained by solution for AL of the equation 


f 0 T = g^ 2 / 202 


( 12 ) 
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As an estimate of f 0> the value used in reference 5 was the fre- 
quency of first symmetrical wing tending, the dominant mode in wing 
buffeting. This frequency was 11-7 cps for run B. The 'probability 
distribution of the peak wing loads of run B, expressed as L/P, iB 
plotted in figure 10, in comparison with the asymptotic expression for 
these loads. The lines represent three different values of fQ. The 
line marked 11. 7 cps corresponds to the assumptions of reference 5> 

14.8 is the value of fQ counted (table IV), and 16.4 cps is the spec- 
trum value (eq.. (4)). Each of the lines is a fair representation of the 
data. The differences between them would be of little practical signif- 
icance in such fields as fatigue, where differences of 2 to 1 in such 
quantities as fatigue life are to be expected. 

For the random response of a linear single-degree-of -freedom oscil- 
lator subject to a Gaussian random input with a white spectrum, the zero 
crossing frequency fg is equal to the undamped natural frequency. The 
difference between 9-5 cps for run A, 11. 7 cps for run B, and the values 
of fg given in table IV may be considered as the effects of the pres- 
ence of other structural modes in the buffet loads (fig. 6) and of 
departures of the color of the unknown buffet input spectrum from white; 
nonetheless, it appears from figure 10 that a single-degree-of -freedom 
approximation would account at least for first-order effects. 

In the present paper, use of the root-mean-square buffet loads over 
short intervals, as in figures 2, 5# and 4, has resulted in efficient 
use of flight-test results. In the analysis of reference 5, only one 
datum point was obtained for each run, since the measure of buffet load 
was AL, one-half of the largest peak-to-peak fluctuation measured in 
the run. The statistical nature of the relationship between the standard 
deviation or root mean square and the peak values makes possible a com- 
parison of the results obtained with the two measures . In 'reference 5 
the analysis of the wing loads in 26 runs led to the relation 

= (155-5 ± 6.4)Pfi (13) 

The slope of the least-squares line through the wing data of figure 4(b) 
is 672 + 25 pounds per unit change in penetration which, taken with a 
value for q. of 292 pounds per square foot, leads to 

0*7 - 59P (1^) 

The 26 maneuvers analyzed in reference 5 averaged 4.6 seconds in 
duration. If fQ is taken as 14.8 cps (table IV ) , equation (12) indi- 
cates that on the average not more than one peak in 4.6 seconds' would 
exceed the standard deviation by a factor of 2.92, or 
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~ 2.92cfy 

« LUff ^ (15) 

The results for run B are therefore within about 30 percent of the results 
obtained for the 26 runs of reference 3 • The agreement would have been 
closer had the values of P used for the correlation of reference 3 been 
the largest penetration rather than the sometimes smaller actual value 
at the time the peat load was obtained, but the results frcm the two 
different measures of buffet intensity can be considered in good agree- 
ment. Whereas the standard deviation mates more efficient use of the 
data and utilizes a measure of greater stability, more computational 
effort is required to obtain it than the simpler peat-to-peat values. 
However, either technique could be used, the standard deviation being 
the preferred measure when flight-test time. must be conserved. 


CONCHJDING REMARKS 


Same statistical properties of the buffet loads measured on the wing 
and tail of a fighter airplane have been analyzed in an attempt to deter- 
mine whether buffeting can be considered a Gaussian random process. 

For a representative stall maneuver it appears that the wing buffeting 
is essentially a Gaussian random process. The loads are normally distri- 
buted, and the probability that a load peak will exceed a given level is 
in agreement with the theoretical results obtained by Rice in "Mathematical 
Analysis of Random Noise. " There is evidence that the tail buffet loads 
are not so normally distributed as the wing loads, and there is also same 
evidence of correlation between angle of attack and tail buffet intensity. 

For a representative pull-up into buffeting in the Bhoek regime, the 
buffet intensity appears to vary linearly with penetration beyond the 
buffet boundary. The loads under maneuvering conditions are therefore 
not stationary and are thus non-Gaussian, but by means of a simple linear 
transformation the buffet loads in maneuvering flight can be treated as 
a Gaussian process. 

The power spectrum of the wing root shear Indicates that buffet 
shear loads are primarily associated with response in the first symmetri- 
cal bending mode. Although other structural modes are present, their 
contribution is small, and the spectrum is essentially characteristic of 
the response of a lightly damped single-degree-of -freedom system to a 
random disturbance. The power spectrum of the tail root shear indicates 
that the tail loads are primarily associated with the fuselage in a tor- 
sion, or "tail-rocking, " mode. 
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In studies of the relationship between buffeting intensity and air- 
plane operating conditions, numerical filtering techniques have been 
found effective for separating the buffet loads from the low-frequency 
load components associated with the maneuvering of the airplane. A use- 
ful measure of buffet intensity was the root -mean-square buffet shear 
averaged over a l/2-second interval. This measure was found to be suffi- 
ciently stable statistically, yet did not encompass so long an interval 
that the effects of changes in flight condition were obscured. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 10, 1956. 
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APPENDIX A 

SOME CHARACTERISTICS OF A GAUSSIAN RANDOM PROCESS 


The nature of some of the concepts available for the description 
of random processes and the simplifications which result when the pro- 
cess is Gaussian are discussed in this appendix. Eknphasis is given to 
some of the practical applications of items which are developed in more 
detail in references k, and- 7. The discussion is restricted to 

what is termed a "stationary 1 random process . For random time series 
such as are dealt with in the present paper, the restriction to stationary 
means that while the quantity of interest (for example, load) varies with 
time, the value of any statistical measure of the quantity (for example, 
its mean square) does not depend upon the time for which the measure is 
determined. Thus time does not enter directly, hut is used only as a 
means of specifying the duration of a time interval. 


Probability Distributions 

A stationary random function of time y(t) may be characterized 
by a series of related probability functions W-^, W 2 , . . . W n where 


W^(y) dy is the probability of finding y between y and y + dy; 


W 2 (yi,y 2 i T ) dy-^ dy 2 is the joint probability of finding a pair of values 
of y in the ranges y 1 to y-^ + dy 1 and y 2 to y 2 + dy 2 , which 
are a time interval r apart from each other; 

Wj(yi, y 2 , y 5 j' r 1 ,T 2 ) dy-j^ dy 2 dy^ is the joint probability of finding a 
triple of values of y in the ranges y-^ to y 1 + dy^, y 2 to 
y 2 + dy 2 , and y^ to y^ + dy^, where dy^ atJ d- dy 2 are the ti m e 
interval apart and dy 2 and dy^ are the time interval t 2 apart. 


The W n probability distributions represent a complete statistical 
description of the process, each distribution of higher order describing 
it in greater detail. The function W-j_(y) is termed the first probabil- 
ity density; it fulfills the requirements that a probability function is 
never negative, a n d that 



dy = 1 


(Al) 
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It also serves "bo define the probability function P(y), which is the 
probability that a given value of y will he exceeded; that is 



From W-^(y) the average value of y can also he found: 


y = f y ^i(y) ay 

—00 


(A 2) 


(A3) 


which (for a stationary random process) is the same as the time average 



A measure of the spread of the values of 
value y is the standard deviation tf, which 

root of the average value of (y - y) 2 : 



y(t) around the average 
is defined as the square 


(A) 


The value of a may also he obtained from W-jXy) .as 


n oo 

a 2 = / (y - y ) 2 W L (y) dy 

^ — OQ 


(A5) 


For many purposes it is convenient to express y in terms of its 
fluctuations about the average value, normalized by the sta n dard devia- 
tion. Such a dimensionless expression is denoted by the standard vari- 
able z, where 


z 



(A 6) 


When r»hgmgrfng variables from y to z, a dual change of scale is 
involved with probability densities, since the fundamental requirement 
of a probability distribution must be observed; that is. 



(A7) 
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and hence 


W;l(z) = W x (y) g = a W^y) (AB) 

From the second probability distribution, W 2 , is obtained the mean 
value of the product y(t) y(t+T) where t is the difference in time 
between the values of y. This mean value 

y x y 2 = JJ y-]y 2 w 2( y i* y 2' T ) ^2 ( A 9 ) 

which is a function of r, is the same (for a stationary process) as the 
time average 


/ T 

y(t) y(t+r) dt (A10) 

T 

It gives a measure of the correlation between values of y separated 
by a time interval r and is termed the correlation function R(t) of 
the random process . From equation (AID ) when t -= 0 it is evident that 

R(0) = y2 (AH) 

For conformity with standard statistical practice where perfect corre- 
lation is denoted by a value of 1 , the correlation function may be 
normalized and expressed as 


p(t) 


[y(t) 

-y] 

[y(t-H-) - y] 


y(t) 

-y] 



(A12) 


Power Spectrum 

The information contained in the W n probability distribution 
functions is also contained in a set of distribution functions which 
describe the frequency content of the time variations of the random 
function y(t) . The simplest of these frequency distribution functions 
is generally called the power spectrum or power spectral density and in 
this paper is denoted by the symbol ©(f) - Under suitable limitations 
(see, for example, ref. 6 , sec. 6 . 7 ), ®(f) can be defined in terms of 

a Fourier transform of the random time function y(t), but frequently 
of more practical interest is the relationship with the correlation 
function R(t) expressed in the Fourier cosine transform pair 
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4(f) = 

poo 

/ R(t) cos 2j£fT dr 
0 

(A13) 

R(t) -j 

noo 

/ 4(f) cos 2rtfT df 

(Al^) 


Eq ua tion (A13) is the "basis for most numerical methods of obtaining 
$>(f ). From equation (A14-) when r = 0 a relation follows between power 
spectral density and the mean square: 


f 4(f) df = R(0) = y 2 
Jn 


(A15) 


r rhnn the power spectrum may be regarded as a frequency analysis of the 
mean square value of the random time function. The function 4(f) when 
applied to a random process is regarded as a continuous function, in 
dis ti nc tion to, say, Fourier coefficients which convey similar informa- 
tion about the amplitudes of harmonic functions. In the event that a 
time function cont ains nonrandom elements such as a mean level (in elec- 
trical terms, a direct ^current component) or periodic components such 
as A k sin (2irf k t - a^), the spectrum will contain discontinuous peaks. 

This information can, however, be expressed in continuous form through 
the convention of the Dirac delta function &(f-f k ) , defined as 


8(f-f k ) =0 (f i ffc) 



> 


(A16) 


Also, 8(f) = 8(-f ) and hence 



(A17) 


With this convention, for a time function 


y(t) = y + A k sin(2irf k t - a^) 


(A18) 
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tlie correlation function is 

A 2 

R(t) = y 2 + cos 2jtf k T (A19) 

and the power spectrum is 

A 2 

®(f) = 2f%(f) + 6(f-f k ) (A20) 

These relationships are frequently of use in analyzing the characteristics 
of experimentally derived spectra. 


The Gaussian Random Process 

A stationary random process is termed Gaussian provided all members 
of the set of multidimensional probability distributions W n are Gaussian 
This requirement means that the first probability distribution Wp(y) 
may be expressed functionally as 

- jfgj e-fr-y) 2 /* 2 (A2i) 

This equation represents a one-dimensional Gaussian distribution, the 
so-called normal distribution of statistics. The second probability 
distribution W^y^yg, - !") is a two-di m ensional Gaussian distribution 

which involves the normalized correlation function p(t): 

y 1 2 +y 2 2 - 2 Ry 1 y 2 

W**) = 2*^1 - pg 6 202(l ' p2) (A22) 

and in general the W n probability distribution is an n-dimensional 
Gaussian distribution. 

The definition of a Gaussian random process in terms of such a set 
of n-dimensional Gaussian distributions, although complete, is hardly a 
definition of practical utility. Its importance stems from the fact that 
for the Gaussian random process all of the distributions W n will depend 

only on a 2 and p(t), two quantities which in turn are derived from the 
power spectrum, as in equations (A13) and (A15). Therefore, the spectral 
density $(f) may, with equation (A2l), be regarded as a basic descrip- 
tion of a Gaussian random process. Other characteristics of Gaussian 
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process follow from the relationship between the power spectrum and 
the W n Gaussian probability distributions. Of considerable practical 
interest in this connection are the relationships between the power 
spectrum, the derivatives of y(t), the frequency of zero crossings, and 
the frequency of occurrence of peaks in the time history. 

Derivatives of a Gaussian random process .- A consequence of the 
Gaussian character of a random time function y(t) is the relationships 
between the derivatives of y(t) and the power spectral, density 4>(f): 


y(t) S =f^ (2tff) 2 4(f) df = op. 2 (A23) 

and 

y(t) 2 = J (2rtf)\(f) df = c y 2 (A2lf) 


It can also be shown that the correlation between y(t) and its 
derivatives is 


and 


y(t)y(t) = 0 = y(t)y(t) 

> 


(A25) 


y(t)y(t) = -y(t) 2 ' 


For some analytic frequency functions ®(f) which are satisfactory 
representations of certain physical processes over the frequency range 
of usual interest, equations such as (A23) or (A2^) would have no meaning, 
since the integrals would not converge. However, it appears that, for 
physical systems parasitic effects will always be present at sufficiently 
high frequencies, and $(f ) will always approach zero rapidly enough so 
that the integrals will converge. 

Peak: values of a Gaussian process .- In part III of reference t a 
number of results are obtained on the distribution of the zeros, of the 
maximums or peaks, and of the envelope of a Gaussian time function. The 
relations for the number of zeros and for the number and probability 
distribution of the peaks are of special interest, because they are 
easily determined from a time-history record and from the power spectrum 
and because they are quantities which enter directly into the study of 
fatigue. Hie number of times per second that the zero axis is crossed 
with positive (or negative) slope is, on the average. 
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f _ i °y 
T o - — — 


2lt Oy 


and the total number of positive (or negative) peaks is 


fp=^^ 
P 2rt 0- 


(A 26) 


(A27) 


Equations (A26) and (A27), together with equations (A15), (A23) , and 
(A24), lead to equations (4) and (5) of the present study. The ratio 
fg/fp is the negative of the coefficient of correlation between y 

and y. The probability Wp(y) dy that a positive peak will fall 

between y and y + dy is given by a peak probability density function 
which, in terms of the standard variable z, is 


where 


Wp(z) 


*L_ -^K 2 

\fSt 





(A28) 


(A29) 

(A30) 


and is the normal probability that a value 

exceeded} that is. 


z/kg will be 



(A31) 


Equation (A28) is essentially equation (3-6-5) of reference 4 
expressed in the notation of the present paper. The frequency ratio 
fg/fp is the parameter which determines the shape, of Wp(z)} equa- 
tion (A28) is plotted in figure 5(b) for several values of the ratio 
f 0 /fp. For the limit of fg/fp = 0, the distribution reduces to the 

normal distribution} for the upper limit of fg/fp = 1, the distribution 
reduces to the Rayleigh distribution. 
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Equation (A28) has heen integrated to obtain the probability Pp(y) 
that, of all the positive (or negative) peaks of a Gaussian random pro- 
cess, a peal: chosen at random will exceed a given value of y = za. The 
resulting expression 


Fp(y) = / W p (z) dz 
J y/cr 

= ^) + |^ 2/2 f- p N(4)] <«2) 

is plotted in figure 5( a )j it is a function of normal probability distrib- 
utions and the ratio f 0 /fp. 

The actual number of positive peaks per second which would exceed 
a given value of y = za is, on the average, E z , where 

N z = f p P p (z) (A33) 

As pointed out in reference k, an asymptotic expression for the 
peak probability density, equation (A28), is 

Wp(z) « f 0 ze“ z2 / 2 (A34) 

an expression which is quite accurate for z > 1 or y > a. The corre- 
sponding expression for the number of positive peaks per second which 
exceed a given value of z is, from equations (A32) and (A33), 

% « f 0 e" z2 / 2 (A35) 

This is a convenient expression which plots as a straight line on semi- 
logarithmic paper when z^ is used as the abscissa. 
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APPENDIX B 
NUMERICAL FILTERING 


The load time histories recorded in the present investigation can 
he considered to he composed of a maneuver component and a "buffet com- 
ponent- The frequency content of the maneuver wan concentrated in the 
vicinity of the frequency of the short-period pitching motion (about 
1 cps) while the principal components of the buffeting were of the order 
of 10 to 12 cps. For separation of the buffet component, same method of 
discrimination between components above and below about 3 to 4 cps was 
desired. Estimation of the maneuver component by means of a manually 
faired mean line and measurement of the buffet component relative to this 
line was tedious and uncert ain . Estimation of the maneuver component by 
calculation of the running mean proved to be a simple numerical procedure, 
but use of the r unning mean is equivalent to use of a low-pass numerical 
filter with a long tail and systematic phase shifts which were undesirable. 
A special high-pass numerical filter was therefore devised which is suit- 
able for use with punched-card equipment, has no phase shifts, and has a 
sharply defined cutoff in the desired frequency range. The filter was, 
in effect, based on the properties of a particular class of low-pass, 
high-pass, and band-pass numerical filters. Since the approach used can 
be extended to the development of other filters with a wide range of 
characteristics, and since it differs in seme respects from the method 
used in reference 9, it is described in seme detail. 


Frequency Characteristics of Numerical Data 

For application of numerical procedures, a continuous function of 
time y(t) is ordinarily reduced to discrete time-history form y^ by 
sampling at the uniform time interval At. This new function y^ is 
thus defined at the particular values of 

t = i At (i = 0, 1, 2, . . .) (Bl) 

but is undefined at intermediate values. There is no loss of information 
in this reduction if y(t) contains no components with frequencies 
greater than a frequency fp, termed the Nyquist frequency or folding 
frequency, where 



1 

2 At 


(B2) 


The folding frequency has the property that in any frequency representa- 
tion of y(t), dete rmin ed from the values of y^, components of 
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frequency f, (2f F ± f), (4f F ± f), and. so forth, cannot he distinguished. 
Thus components of such frequencies greater than fp appear in the fre- 
quency representation at the frequency f in the range ‘0 5? f < f^,, and 

the whole range of frequencies in y(t) has been folded into the range 
0 to fp. This folding property follows from the relations 

sin 2rt(2f-p ± f)t = sin(2rti ± 2rtft) = ±sin 2rtft 

(i=0, 1, 2, . . .) 


Binomial-Coefficient Filters 
The binomial coefficients as displayed in the array 


k 

n 

0 

1 

2 

3 

h 

5 

6 7 8 ... 

1 

1 

1 






2 

1 

2 

1 

• 




3 

1 

3 

3 

1 




k 

1 

h 

6 

4 

1 



5 

1 

5 

10 

10 

5 

1 


6 

1 

6 

15 

20 

15 

6 

1 

7 

1 

7 

21 

35 

35 

21 

7 1 

8 

1 

8 

28 

56 

70 

56 

28 8 1 

• 

n 

• 

1 

T 1 

• 

r 

i(n- 

• 

l)(n ■ 

- 2). 

• • • 
. .(n - k + l) 



XX 


k(k - 

■ 1 )(k 

- 2) 

...(3)(2)(1) 

and their 

sums 

2 31 form the basis 

for 

a number of integrating 1 


n - 1 


formulas. The smoothing action on numerical data is aiHn to the action 
of filters on electrical signals, as may be seen, for example, by consid- 
ering the action of the formula based on (^ji 
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y(t) = -^|y(t-^t) + ?y(t) + y(t-hAt )] 


(B3) 


when y(t) is a sine wave of frequency f, that is 

y(t) = sin 2rtft 

or, if t = i At, 

. itfi 

y± = sin 3— 
r F 

Application of the relation sin(a + h) = sin a cos b ± cos a sin b 
to equation (B3) gives 


y(t) * 


_ 1 + cos(2jCf At) 


sin 2rtft 


or 


y ± = 


1 + cos 
% 


y * = ( oos2 1 |) y i 


Thus y(t) (or y^) is a sinusoid of the same frequency as y(t) (or 

yi) but the amplitude has been modified by the term cos 2 % which 

2 f p 


thus represents the filter factor of the 


(I) 


binomial-coefficient filter. 


For data-processing purposes, the operation represented by equa- 
tion (B3) is more conveniently expressed in matrix form as 


* 1-1 


*i = J L 1 2 L J ^ yi 


*i+l 


(B4) 


or, in general, for the nth-order binomial coefficients, 

yi = M*)J W 

and the filter factor of the associated filter is easily shown to be 

(B5) 


r n (r) 
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Application of a second-order filter to obtain y^, followed by a repeti- 
tion of tbis action on y^ to give y^, is equivalent to a single appli- 
cation of a fourth-order filter, as may be seen from the’ matrix repre- 
sentation of this double filtering process: 





’* 1-1 
- yi ► 
^i+l 



0 12 10 
0 0 12 1 


^i-2 

. yi-i 

- y ± ► 

y i+l 

y i+ 2 


== I L 1 ^ 6 ^ iJ {yi} 


Successive operations of this type are commutative, use of an nth-order 
filter followed by an mth-order filter being equivalent to use of an 
(n+m) -order filter, and the filter factors are also commutative; that is. 


, »n-fm 

Y (n«)M - = (“ | 4 ) (B6) 


The binomial filters are all low-pass filters, the filter factor 
being unity at f = 0 and zero at the folding frequency. The filters 
determined by the coefficients of even order have the special properties 
that they are nonnegative and have no phase shifts; thus, the time 
relationships of components of different frequencies are preserved. The 
low-pass -filter characteristic Y n (f) is plotted in figure ll(a) for 

the range 0 S S- ^ 1 for several even values of n. The higher the 
r F 

order of the filter, the sharper the low-frequency cutoff. As measured 
by the usual criterion of band width at the ha lf -power point, Y(f) = O. 7 O 7 , 
the second-order and 28 th-order filters would have band widths of 0 . 36 f -jr 

and O.lOfy, respectively. 


Modified Binomial Filters 

By a simple operation on the coefficients, each even-order binomial 
filter can be related to a high-pass filter which has a filter factor 
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of 0 at f = 0 and 1 at the folding frequency. The time history z^ 

computed with the high-pass numerical filter is related to the time history 
y^ by the relation 


z i = yi - y± 


(B7) 


that is, the double operation - application of a low-pass filter followed 
by substraction from the original time history - is equivalent to a single 
operation with the related high-pass filter. 

Since, from equations (B4-) and (B7), 


y i-(n/2) 


z i = yi 





y i+(n/2>J 


the coefficients for finding in one operation are contained in the 

relation 


y i-(a/2)] 



y i+(n/2) 

. ■* 

(B8) 

That is, the coefficients for the nth-order high-pass filter are the 
negatives of the coefficients for the low-pass filter, except for the 

central coefficient k = S, -which is equal to 2* 1 minus the central 

binomial coefficient ^ n ^>) • Htus the sum of all the coefficients of the 

high-pass filter is zero. In the following arrays the left-hand set 
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represents the low-pass filters and the right-hand set the related high- 
pass filters. The normalizing factor for each set of coefficients is 
the sim of the binomial coefficients, 2? 1 . 


n 

2 11 

Low Pass 

High Pass 


2 

h 

1 

2 1 

-1 

2 -1 


k 

l 6 1 

k 

6 k 1 -1 

-4 

10 -4 

-1 

6 

6ii- i 6 

15 

20 15 6 1 -l -6 • 

-15 

14 -15 

-6 -1 

The high-pass-filter characteristic Z^f ) _ for 

the 

nth-order 

filter is 




Zn(f) = 1 - Y n (f) 







= 1 - cos n | 

2 r p 



(B9) 


which, for n = 2 , is sin 2 ~ The filter characteristic Z^f) is 

d. Xji 

plotted in figure ll(b) for comparison with Y n (f) (fig. 11 (a)). As 
the order increases (and thus the length of time over which the filter 
"remembers" increases), the suppression of high-frequency components by 
the low-pass filter becomes more effective, . and the filter becomes 
sharper j similarly, the pass band of the high-pass filter becomes wider. 

Operations with low- and high-pass filters are commutative, and give 
band-pass filters. For example, combination of the second-order low-pass 
filter jjsjl 2 ij and the high-pass operator [-1 2 -ij is equivalent 

to the filter ^|_-1 0 2 0 -lj, which has a characteristic 

F(f ) = lff 2 (f)Z2(f) 

= sin 2 0 - (BIO) 

f F 

Similarly, the operator i[-l 0 lj (based on the first-order binomial 
coefficients) has a filter amplitude characteristic of 


F(f) = sin 

r F 


(Bll) 
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The "band-pass filters can also "be regarded as modified "binomial 
filters for which the folding frequency has "been redefined. For example, 

the operation j^|_“l 0 2 0 -ij can "be regarded as the high-pass fil- 
ter of order 2 operating on the points y i-2 , Y±> and- y±+2> tlius > 111 

effect, the time interval At "between data points is doubled, and the 

f-m 1 

effective folding frequency is the value -g- = Thus if 

jj=|j-l 2 -ij is represented "by a filter factor 

Z 2 (f) = sin 2 £ f- 

2 f p 

the operator j^[-l 020 -lj represents the filter factor 

F(f) = sin 2 | (B12) 

That is, 

F(f ) = sin 2 (B13) 

r F 

which is the same as equation (BIO). 


Separation of Maneuvering and Buffet Loads 

For the purposes of the present study a pair of low- and high-pass 
filters was developed. These filters, like the 28th-order binomial fil- 
ters, operate on 29 successive values of the measured time history y^, 
but differ from the binomial-type filters in both sharpness of discrimina- 
tion at low frequencies and in suitability for routine computations with 
punched-card equipment. The filters are more easily described in t erms 
of the low-pass member of the pair, which was constructed from the low-pass 
fourth-order binomial operator 

^|_1 k 6 k lj (Bl&a) 

and the two related modified binomial operators with band-pass 
characteristics 


J^[l 0 lj. 060 lj. 0 lJ 


(Bllj-b) 



32 


NACA IN 3733 


J-|l 000400060004000 

1 6 L 


Ij 


(Bl 4 c) 


These operators, respectively, have the filter characteristics 


%( f ) ■ I £ 


Y^(2f ) = cos 4 | ft 
?4(4f) = cos 1 *- * 

F 

as shown in figure 12 (a). The three operators (Blit-) can he used in 
succession with three passes of the simpler punched-card calculators, 
or with more elaborate equipment the filtering can he accomplished with 
the combined symmetrical operator which acts on 29 adjacent points: 


r— -li it 10 20 35 56 84 120 161 204 2*4.6 284 315 336 344 336 . . . 1| 
4096 L 

(B 15 ) 

In either case, the filter factor F(f ) is given hy the product of the 
filter factors of the three operators as 



which is illustrated in figure 12 (h) . As measured hy the criterion of 
hand width at the half -power point, F(f) = O.7O7, this filter has a 
hand width of 0.0575f F , or as used in the present studies, where 

f F = 50 cps, the hand width is 2.9 cps. It is thus sharper than the 

28th-order binomial filter (fig. 12(a)) which operates on the same num- 
ber of points. Four-tenths of one percent of the total transmission 
lies in three minor lobes located between the zeros at X_ = 0.25, O.5O, 

f F 

O.75, and 1.00, and shown to an expanded scale in figure 12(h), hut these 
lobes are all smaller than 0.0028 in amplitude. 

In order to obt ain the buffet load, the time-history values of the 
man euver load obtained through use of the operator (B 15 ) were substracted 
from the measured load time history. It could have been obtained directly 
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■by use of the high-pass filter equivalent to the operator (B15) hut, 
since the central coefficient of the filter would have then contained 
not three hut four digits (4,096 - 344 = 3,752), the "buffet load was 
obtained in two steps in order to preserve the computational efficiency 
associated with the three-digit coefficients of (B15) . 
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TABLE I.- GROUND VIBRATION FREQUENCY OF AIRPLANE STRUCTURAL MODES 



Frequency, 

cps, for • 


Basic 

airplane 

Modified 

airplane 

Wing modes: 

First symmetric tending 

First asymmetric tending ' 

First torsion 

Second symmetric tending 

11.7 

38.0 

9-3 

18.1 

34.5 

52.0 

Horizontal statilizer modes: 

First symmetric tending 

First asymmetric tending 

First torsion 


25.0 

36.0 
70.0 

Fuselage modes: 

Torsion 

Side tending 

Vertical tending 

9.8 

12.5 

9.-8 

12.5 

ih.9 




Bun A 



Frequently for - 
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TABLE IH.- FREQUENCY DISTRIBUTION OF TRANSFORMED 
LOADS AND LOAD PEAKS FOR RUN B 




Frequency for - 

Range of 

TK 

L/P, 

Wing 

Tail 

J-D 










Loads 

Positive 

Negative 

TnnrlR 

Positive 

Negative 




peats 

peats 


peats 

peats 

-2,000 to • 

- 1,800 

1 


1 




- 1,800 to 

- 1,600 

3 


3 




-i,6oo to 

- 1,400 

3 


2 




-1,400 to 

-1,200 

5 


2 

2 


2 

-1,200 to -1,000 

ID 


6 

4 


4 

-1,000 to 

-8oo 

11 


8 

10 


6-. 

-800 to 

-6oo 

21 


8 

22 


15 

-600 to 

-4oo 

23 


13 

30 

1 

8 

-400 to 

-200 

38 

2 

7 

42- 

4 

16 

-200 to 

0 

36 

6 

3 

40 

3 

9 

0 to 

200 

40 

6 

2 

47 

7 

6 

200 to 

4oo 

4i 

8 

2 

50 

14 

5 

400 to 

6oo. 

27 

6 

2 

. 35 

23 

l 

600 to 

8oo 

21 

12 


21 

7 

1 

800 to 

1,000 

20 

13 


7 

5 


1,000 to 

1,200 

10 

6 


6 

5 


1,200 to 

1,400 

5 

3 


2 

1 


1,400 to 

1,600 

3 

3 


1 

1 


1,600 to 

1,800 

1 

1 


0 

0 


1,800 to 

2,000 

0 

0 


1 

1 


2,000 to 

2,200 

1 

1 





N 

320 

67 

67 

320 

73 

72 

Mean 


25 



25 



cr 


653 



512 



X 2 


19.8 

4.7 

4.8 

5-5 

2 

x .05 


22.4 

11.1 

14.1 

11.1 
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TABLE IV.- VALUES OP FREQUENCY PARAMETERS' OBTAINED 
FROM POWER SPECTRA AND FROM TIME HISTORIES 


Source 

Wing 

Tail 

f 0> 

cps 

f pj 

cps 

f o/ f p 

fo> 

cps 

f P> 

cps 

f o/ f p 

Run A: 

Spectrum 
Time history 

18.1 

16.3 

00 CM 

0.52 

.62 

15.5 

15.7 

25.4 

20.0 

0.61 

.78 

Run B: 

Spectrum 
Time history 

16 .4 
14.8 

27.3 

21.8 

.60 

.68 

15.5 

16.1 

28.1 

22.8 

•55 

.70 

Run B (transformed): 
Time history 

14.8 

21.0 

•77 

16.1 

22.7 

•71 




r 


a 


(a) Run A. 


Figure 1.- Maneuvers analyzed in present investigation. 
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(a) Variation with Cjj. 

Figure Variation of wing and tail "buffet 



(b) Variation with penetration P. 
load with Cjj and P in shock regime. Run B. 
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(a) Probability Pp(z) that a peak will exceed a given value. 

Figure 5*- Peak probability functions of a Gaussian random process for 

various values of the ratio fg/fp. 
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(b) Peak probability density Wp(z). 
Figure 5»- Concluded. 






of buffet load 
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Figure 6. 
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Figure 8.- Probability that a value of transformed load L/p will be 

exceeded . 




Probability of exceeding, 
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h9 



Standard variable, z 


Figure 9*- Probability that a peak value of buffet load will exceed a 

value equal to z a. 








Filter factor, Y n (f) 


W 


f 


| 

I 

J 


) 



’ |W ** »w IV JL|V V I" • * ; — — * — 

Frequency ratio, f/fp Frequency ratio, f/i^ 

(a) Low-pass filter. ("b) High-pass filter. 

Figure 11.- Filter factor of 'binomial filter of order n and associated higgi-pass filter. 


I 


cac hi vom 



- LuntUy Field, V*. 


VJ 1 

ro 




(a) Elementary filters. 


(b) Low-pass filter vised. 


™j T ft iVp 

J _H£lir« IC*" V/UilQblLIUbXUll UJ. 


4>Via «v.mia<w4 nel •P^l + ov -PriT* flOfnOm+. l TICT +.n+Mj. *1 OAfi TT1A Afluremfiirt B into 

MiQ » ttUMB ii -«*. 1 -*. . I .j . u\«d. *W* MV jy w ** w,>, ***0 w »» — — . 

maneuvering and buffet components. 
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